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Irreversible aggregation processes involving reactive and frozen clusters are investigated using the 
rate equation approach. In aggregation events, two clusters join irreversibly to form a larger cluster, 
and additionally, reactive clusters may spontaneously freeze. Frozen clusters do not participate in 
merger events. Generally, freezing controls the nature of the aggregation process, as demonstrated 
by the final distribution of frozen clusters. The cluster mass distribution has a power-law tail, 
Fk ~ k when the freezing process is sufficiently slow. Different exponents, 7 = 1, 3 are found 
for the constant and the product aggregation rates, respectively. For the latter case, the standard 
polymerization model, either no gels, or a single gel, or even multiple gels may be produced. 
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I. INTRODUCTION 

Aggregation is a fundamental irreversible process in 
which physical objects merge irreversibly to form larger 
objects. Aggregation has numerous applications ranging 
from astronomy where planetary systems form via colli- 
sions of planetesimals, to atmospheric science 0, j to 
chemical physics where polymeric chains chemically bond 
and forrnpolymeric networks or gels 0,0,13 j to computer 
science laHj^. 

The standard framework for modeling aggregation is 
as follows. Initially, the system consists of a large num- 
ber of identical molecular units ( "monomers" ) . A clus- 
ter ("polymer") is composed of an integer number of 
monomers, termed the cluster mass. In each aggregation 
event, a pair of clusters merges, thereby forming a larger 
cluster whose mass equals the sum of the two original 
masses. 

When the number of aggregation events is unlimited, 
the system condenses into a single cluster. However, 
in most practical applications, other processes intervene 
well before that, and as a result, the final state has mul- 
tiple clusters, rather than a single condensate. For exam- 
ple, fragmentation of large clusters into smaller clusters 
is one mechanism that may counterbalance aggregation 
and prevent condensation. 

In this study, we focus on another control mechanism: 
freezing. We use the generic term freezing to describe 
situations where there are two types of clusters: reac- 
tive clusters that participate in aggregation events and 
passive clusters that do not participate in aggregation 
events. 

In particular, we consider the case where reactive clus- 
ters have a finite lifetime. In our model, reactive clus- 
ters spontaneously turn into frozen clusters. Spontaneous 
freezing can occur via various mechanisms. For instance, 
the environment may contain "traps" that absorb the 
diffusing polymers. In this situation, the reactive clus- 
ters are the free (mobile) polymers and the frozen clus- 
ters are the polymers adsorbed to the trap surface. An- 



other example is a system of linear polymers with reactive 
end monomers. In a merging event, two different chains 
chemically bond via the end monomers, while in a freez- 
ing event, the two end monomers of the same chain bond 
to form a ring. Rings can no longer participate in aggre- 
gation events. Thus, in this case the linear polymers are 
reactive and the ring polymers are frozen. 

In aggregation with freezing, it is natural to consider 
the initial condition where there are reactive clusters 
only. Of course, the system ends with frozen clusters 
only. Of special interest is the final mass distribution of 
frozen clusters. In this study, we address the two classical 
aggregation rates. 

First, we study the simplest aggregation process where 
the aggregation rate is independent of the cluster mass. 
We find that the mass distribution of reactive clusters de- 
cays exponentially with the cluster mass. In general, the 
mass density of frozen clusters also decays exponentially 
with the cluster mass. However, when the freezing rate 
is very small, there is a power-law behavior, Fk ~ k~^, 
over a substantial range of masses. 

Second, we consider the case where the aggregation 
rate is proportional to the product of the masses, a pro- 
cess that is widely used to model polymerization and 
gelation. We find that when the freezing rate exceeds a 
certain threshold, no gels form, while when the freezing 
rate is below this threshold, at least one gel forms. In- 
terestingly, the number of gels produced fiuctuates from 
realization to realization. For supercritical freezing rates, 
the mass distribution of frozen clusters decays exponen- 
tially, while below this threshold, it decays algebraically, 
Ffc ~ fc-3 i. 



II. THE MASTER EQUATIONS 

We analyze the stochastic process of aggregation with 
freezing using the rate equation approach. Let us first 
consider the evolution of reactive clusters. In aggrega- 
tion processes, two reactive clusters of masses i and j 
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merge to form a larger reactive cluster of mass i+j. The 
aggregation rate K{i,j) is a function of the two cluster 
masses. The freezing process is random: reactive clus- 
ters may spontaneously freeze with a constant rate. This 
freezing rate fk may be mass-dependent. Therefore, the 
mass distribution Rk {t) of reactive clusters of mass k at 
time t evolves according to the generalized Smoluchowki 
equation 



dRk 
dt 
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The first two terms account for gain and loss of clusters 
of mass k, and the last term accounts for loss due to 
freezing. This master equation assumes perfect mixing 
as the probability of finding two clusters at the same 
position is a product of the probabilities of finding each 
of the clusters independently at the same position. We 
restrict our attention to the natural case of monodisperse 
initial condition, i?fc(0) = Sk,Q- 

The mass distribution of frozen clusters Fk{t) is cou- 
pled to the mass distribution of reactive clusters accord- 
ing to the rate equation 



dFk 
dt 



fkRk- 
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It is simple to check that the total mass density 
J2k f^i^k + Fk) = 1 is conserved by the evolution equa- 
tions O"©- Initially, there are no frozen clusters so 
Ffe(O) = 0. Eventually, all clusters become frozen, so 
the final mass distribution Fk{oo) of frozen clusters is of 
particular interest. 

The master equations CQ"© are sets of infinitely 
many coupled nonlinear differential equations, and they 
are generally unsolvable. Even in the absence of freez- 
ing, these equations are solvable only for special aggre- 
gation rates P, 0, 0| . The three classical solvable cases 
are the constant rate K{i,j) = const., the sum rate 
K{i,j) = i+j^ and the product rate K{i,j) — ij. These 
cases represent natural aggregation processes. Mass in- 
dependent aggregation rates correspond to an aggrega- 
tion process where two clusters are chosen randomly to 
merge. Aggregation rates proportional to the product of 
the two cluster masses correspond to polymerization pro- 
cesses where two monomers are picked randomly to form 
a chemical bond; consequently, their respective clusters 
are merged. The sum rate is a hybrid between the two 
as it is an aggregation process where a randomly chosen 
monomer bonds with a randomly chosen polymer. In 
this study, we focus on the two most widely used cases 
of constant and product aggregation rates. 



equations with inhomogeneous densities and add diffu- 
sion terms. Then, one should study these equations in 
the trap-free region subject to the absorbing boundary 
conditions imposed by the traps. This approach is not 
practical and the reaction-rate approach provides a pow- 
erful alternative 0, ITl. IT^ Il5j | . The reaction-rate 
approach is roughly speaking an effective-medium the- 
ory that ignores the complicated influence of each trap 
on the diffusion of particles and instead, represents this 
influence by averages. The reaction-rate approach was 
used by Smoluchowski to compute the aggregation rate 
K(i,j) for Brownian particles. Assuming that merging 
happens immediately upon collision, and that particles 
are spherical and have radii Ri and Rj and diffusion co- 
efficients Di and Dj, Smoluchowski obtained 

K{i,j) = 47r(A + Dj){R, + R,). (3) 

Stokes' law shows that the diffusion coefficient of a Brow- 
nian particle is inversely proportional to its radius, Dk ^ 
1/Rk ~ k^^^^, and therefore the Brownian kernel be- 
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Here, we ignored the overall multiplicative factor as it is 
irrelevant for the current discussion. 

The master equations with this complicated Brownian 
kernel have not been solved even in the case of pure aggre- 
gation. To simplify the analysis, Smoluchowski suggested 
to replace the Brownian kernel Q by the constant kernel. 
These two kernels have one common feature — they both 
are invariant under the dilatation K(ai,aj) = K{i,j). 
Therefore, one expects that both kernels lead to similar 
behaviors, and to a certain extent, i.e., as far as overall 
scaling properties are concerned, this approximation is 
sensible p^. 

A straightforward extension of Eq. ^ gives the freez- 
ing rate 



/fe =4^n iDk + D){Rk 



(5) 



where n is the density of traps that are assumed to be 
spheres of radius a and diffusion coefficient D. The clus- 
ters are usually polymers whose molecular weight is small 
compared to the size of the traps; hence i?^ <C a and 
Dk ^ D. Therefore fk = AiranDk, yielding the mass 
dependence 



(6) 



Thus, a constant aggregation rate together with spon- 
taneous freezing approximate aggregation of Brownian 
particles in the presence of traps. We stress that the 
use of a constant aggregation rate instead of Q is an 
approximation. 



III. CONSTANT AGGREGATION RATE 

First, we discuss how the constant aggregation rate re- 
lates to polymerization in the presence of traps. To treat 
this problem formally, one should write down the master 



A. Constant freezing rates 

Since the constant aggregation rate merely sets the 
overall time scale, we may conveniently set its value 
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K{i,j) = 2 without loss of generality. Let us first con- 
sider constant freezing rates, fk = a. The master equa- 
tion (QJ becomes 



dRk 
dt 



= ^ R,Rj - {2R + a)Rk 



(7) 



i+j=k 



Here, we used the total density of reactive clusters, 
R = J2k^k. In general, for mass-independent freezing 
rates, it is possible to eliminate the freezing term from 
the master equation by transforming the mass distribu- 



tion Rl. = Cj 



k e 



and introducing the time variable 



1 



a 
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The time variable r grows from to a"^ as the physical 
time increases from to oo. With these transformations, 
the governing equations for the densities Ck reduce to the 
pure aggregation case dCk/dr = J2i+j=k ~ 2C Cfe 
with the total density C = J2k^k- We briefly re- 
call how to solve these equations. The total density 
obeys dC/dr = — and subject to the initial condi- 
tion C{0) = 1, the total density is C(t) = + Let 
us now introduce the exponential ansatz C(r) = Aa''~^ 
with ^(0) = 1 and a(0) = 1 to satisfy the initial con- 
ditions. Substituting this ansatz into the master equa- 
tion and equating mass- independent and mass-dependent 
terms separately, yields dA/dr — —2{1 + t)~^A and 
therefore, A = {1 + t)~^, and da/dr = A leading to 
a = t/{1 + t). The well-known solution for the pure 
aggregation case is therefore 



Ck{T) 



^k-1 



(9) 



Thus, the mass distribution of reactive clusters reads 



Rk{T) = (1-ar) 



(1 + ^) 



fc+i 



(10) 



The exponential mass dependence is as in the pure aggre- 
gation case. Also, the total density of reactive clusters is 
i? = (1 — qt)/(1 + r) and as expected, the reactive clus- 
ters do eventually deplete R{t = cx)) — R{t — 1/a) — 0. 

The mass distribution of frozen clusters is found 
by integrating the equation dFk /dr — aCk with 
respect to time. Substituting lO, and using 
dr/dt = (1 — ar) — e^"*, the integration is immediate 
and 



Fkir) 



1 



(11) 



We see that in addition to the dominant exponential be- 
havior, there is an additional algebraic prefactor. The 
total density of frozen clusters F = Y^j-Fk is found by 
summation, F{t) = Q:ln(l -|- r) and in particular, the 
final density of frozen clusters is F{oo) = F{t = oo) = 



aln(l -I- l/a). Also, the final mass distribution of frozen 
clusters is 



1 



fc V 1 + a 



(12) 



In general, the mass distribution decays exponentially, 
but there is a algebraic correction. 



B. Slow freezing 

The most interesting behavior occurs in the slow freez- 
ing limit: as a — > 0, the final mass distribution becomes 
algebraic 



i^fc(oo) 



(13) 



This power law holds over a substantial mass range, 
k <^ . Beyond this scale, the tail is exponential, 
Ffe(oo) ^ak-^e-"^. 

The results in the slow freezing limit can be alterna- 
tively obtained using perturbation theory. Indeed, the 
modified time variable coincides with the original time 
variable, r — ^ i as a — > 0, and the pure aggregation re- 
sults are recovered. In other words, the freezing loss term 
—fkRk can be neglected in the master equation (Q). Us- 
ing this perturbation approach we address two related 
problems: general freezing rates and aggregation in low- 
dimensional systems. 

Let us consider general mass dependent freezing rates 
fk- Dropping the loss rate from the master equation, 
the reactive cluster density is as in the pure aggregation 
case Rk = i''~^(l -I- t)^'^^^, given by Eq. The mass 
distribution of frozen clusters is obtained by integrating 
Eq. 131 



t 



1 + t 



We see that the algebraic prefactor k 
fore, the final mass distribution is 



(14) 



is generic. There- 
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This behavior applies for masses below some threshold 
mass fc*, while the mass distribution sharply vanishes 
above the threshold. The threshold mass is estimated 
from mass conservation: 



l = J2kFk{oo)^J2f'^ 



(16) 
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k=l 



As argued above, for Brownian coagulation in the pres- 
ence of traps, fk — (3k~^/^. For slow freezing, P 1, 
we conclude that the final mass distribution is algebraic, 
Fk{oo) ~ /3/c~^/^, below the threshold mass k^ ~ 

The rate equation approach neglects spatial correla- 
tions as the probability of finding two clusters at the 



4 



same position is represented by the product of the prob- 
abiUties of finding each cluster separately at the same po- 
sition. This mean-field approximation is valid only when 
the spatial dimension exceeds the critical dimension dc 
[l4.. .15] . ft is therefore interesting to study the behavior 
below the critical dimension. 

We address here the Point Cluster Model (PCM) where 
the radii and the diffusion coefhcients are both mass- 
independent. In this case dc = 2. The lattice PCM is 
defined as follows: clusters occupy single lattice sites and 
hop to adjacent sites with rate D; if a reactive cluster 
hops onto a site occupied by another reactive cluster, 
both clusters merge. We assume that frozen clusters do 
not affect reactive clusters. The PCM without freezing 
can be solved exactly in one dimension. When all lattice 
sites are initially occupied by monomers, the density of 
reactive clusters of mass k is [la IitII 



Rk{t) 



-4Dt 



Ik-i{4Dt) ~ h+i{4Dt)] (17) 



where /„ is the modified Bessel function of order n. 
Here, we implicitly considered the slow freezing limit. 
The density of frozen clusters is found from dFk/dt = 
fkRk, that is of course always valid. Using the identity 
dxe~^ [Ik-i{x) — Ik+i{x)] — 2, the final distribution 
of frozen clusters is 



(18) 



Remarkably, the very same answer (|18|l is also found for 
the continuous version of the PCM. The mass distribu- 
tion H18|) holds up to a certain threshold mass. For ex- 
ample, for the constant freezing rate fk = a 1 the 
threshold mass is fc* ~ ^/D/a. 



IV. PRODUCT AGGREGATION RATE 

The product aggregation rate K{i,j) = ij is equiva- 
lent to the Flory-Stockmayer gelation model 0, 0, E| ■ 
In this model, any two monomers may form a chemical 
bond and when this happens the two respective polymers 
become one. Thus, the aggregation rate equals the prod- 
uct of the cluster masses. In this polymerization process, 
a polymer network ("gel") emerges in a finite time, and 
it is giant in the sense that it contains a finite fraction 
of the monomers in the system. Eventually it grows to 
engulfs the entire system. This gelation model is also the 
simplest mean- field model of percolation [T^ f20l| . 

As in the previous section, we analyze in detail mass- 
indcpcndcnt freezing rates, fk = ct, for which the master 
equation ^ becomes 



dRk 
IT 



E 
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ijRiRj 



ikRi. — aRk 



(19) 



where m is the total mass density of reactive clusters. If 
all clusters are finite in size then in = Mi — J2k>i 
where M„ = J^k ^"^fc the general nth moment of the 



distribution. Again, we consider the monodisperse initial 
conditions Rkit) — Sk,i and Fk{0) = 0. 

Low order moments of the mass distribution obey 
closed equations and thus, provide a useful probe of the 
aggregation dynamics. The total mass density of reactive 
clusters satisfies dm/ dt = —am refiecting the loss due to 
freezing, and therefore 



m{t) — m(0)e 



(20) 



The total mass density decays exponentially with the 
physical time, or equivalently, linearly with the modified 
time, m(T) = 1 — ar. Furthermore, the second moment 
includes in addition to the linear loss term, a nonlin- 
ear term that accounts for changes due to aggregation, 
dM2/dt = M| — aM2. Solving this equation with arbi- 
trary initial condition yields 
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1 



(21) 



Divergence of the second moment signals the emergence 
of a gel in a finite time, i.e., the occurrence of the gela- 
tion phase transition. Fixing the freezing rate, the initial 
conditions govern whether gelation does or does not oc- 
cur: gelation occurs when the initial mass is sufficiently 
large, M2(0) > a, but otherwise there is no gelation. 
Conversely, fixing the initial conditions, gelation occurs 
only for slow enough freezing. For the monodisperse ini- 
tial conditions, the critical freezing rate is etc = 1. When 
gelation does occur, the gelation time is 



t„ = --In ( 1 - 



Af2(0) 



(22) 



The gelation point separates two phases. Prior to the 
gelation time, the system contains only finite clusters 
that undergo cluster-cluster aggregation. Past the gela- 
tion point, the gel grows via cluster-gel aggregation. We 
term these two the coagulation phase and the gelation 
phase, respectively. The above expressions for the first 
two moments are valid for the coagulation phase only. 

The mass distribution of reactive clusters is found 
again by transforming the mass distribution Rk = 
e~°''^Ck and the time variable ©• With these trans- 
formations, the problem reduces to the no-freezing case, 
dCk/dt = \ Y.i+j=kij'^iC!3 - kMCk with TO = Me-"*. 
Using the variable u{t) — JJ" M{T')dT' and the transfor- 
mation Ck — Gfer'^^^e^'^", the master equation reduces 
to a recursion equation for the time-independent coeffi- 
cients Gfe: (fc - l)Gk = U2i+j=k'^j'^i^j- This equa- 
tion is solved using the generating function technique. 
The so-called "tree-function" G{z) = J^k kGke^^ satis- 
fies dG/dz — G/(l — G) and the solution of this dif- 
ferential equation obeys Ge^'^ — e^. The coefficients 
Gfc = k ^~'^ l k\ are found using the Lagrange inversion 
formula [2l|. Hence, the mass distribution of reactive 
clusters is 



Rk{T) 



„k-2 



k\ 



(1 - ar) 



k— 1 ^ 



(23) 
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The corresponding generating function ^{z) = 
X^feli kRj-e^^ can be expressed in terms of the tree- 
function 

3?(z) = r"i(l- aT)G(z + hi t-m). (24) 

Exphcitly, the tree function is G{z) — X]/c>i ^^W ^''^ ■ 

The mass distribution (|23|) is only a formal solution 
because the total mass density m and hence the variable 
u are yet to be determined. Prior to gelation, the solution 
can be obtained in an explicit form because the various 
variables are known. From the first moment (|20|) then 
M = 1 and therefore 

u = T (25) 

for t < tg. In this case, the mass distribution decays ex- 
ponentially at large masses and the typical cluster mass 
is finite. When a > 1, there is no gelation transition, 
and this behavior characterizes the mass distribution at 
all times. Otherwise, when gelation does occur, the gela- 
tion time (|22|l is simply Tg — 1. The gelation point is 
marked by an algebraic divergence of the mass distribu- 
tion, Rk{tg) ~ (1 - a)A:~^/^, for large k. 

Using the explicit expression for Rk prior to gelation, 
we can calculate the mass distribution of the frozen clus- 
ters produced up to that point. Substituting l|25(l into the 
formal solution H23|) and integrating dFi^/dt = aRk over 
time using dr/dt — e^"* = (1 — ar) yields the distribu- 
tion Fk{tg) of frozen clusters produced prior to gelation 
[tg = oo for a > 1) 



p., ^ Jl^7(fc,A;) a<l. 



(26) 



where 7(n, x) = dy y"~^ is the incomplete gamma 
function. When a > 1, this quantity equals the final 
distribution of frozen clusters, Fk{oo) = Fk{tg). At large 
masses, the behavior is as follows 



Ffc(oo) 



2 ^ 



a = 1, 
a > 1, 



(27) 



where A = {277)"^^^ a'^ /{a - 1) and B = -I- In a - 1. 
These asymptotic results were obtained using the steep- 
est descent method. Quantitatively, the mass distribu- 
tion differs from that obtained for constant aggregation 
rates. However, qualitatively, there is a similarity: there 
is exponential decay above a critical freezing rate and 
algebraic decay at and below this critical freezing rate. 
For the constant aggregation rate, the critical freezing 
rate vanishes, but for the product aggregation rate the 
critical rate is finite. 

At the gelation transition a giant cluster that contains 
a fraction of the mass in the system emerges. Past the 
gelation point, two aggregation processes occur in paral- 
lel: in addition to cluster-cluster aggregation, the giant 
cluster grows by swallowing finite clusters. Now, reactive 
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FIG. 1: The total mass, the sol mass, and the gel mass versus 
the modified time r for ol = 1/2. 



clusters consist of finite clusters (the "sol") with mass 
s = Ml, and the gel with mass g. The total mass density 
is m = s + g. These three masses are coupled via the 
evolution equations 



dm 
It 
ds 
di 



s{m — s) 
1 — sTe"^ 



(28a) 
(28b) 



The initial conditions are m{tg) ~ s{tg) = 1 — a. The 
first equation reflects that as long as the gel remains re- 
active, mass may be converted from the reactive state 
to the frozen state via freezing of finite clusters. The 
second equation follows from ds/dt — —aMi — g M2, ob- 
tained by summing (|19|l . The second moment is written 
explicitly, M2 = '5l^{z = 0) = s/(l - sre"*), using the 
aforementioned identity G'{z) = G/(l — G). The initial 
conditions are rn(tg) = s(tg) = 1 — a. Once these equa- 
tions are solved, the solution (^5)1 becomes explicit. We 
analyze this equation using perturbation theory in the 
limits a I 1 and a J. as detailed in Appendix A. For 
general freezing rates, we solve these equations numeri- 
cally (Fig. Pi. 

In addition to the freezing of the finite clusters, the gel 
itself may freeze. One way to characterize the gel by its 
maximal possible size gmax = Unit _>oo 3(0- The limiting 
behaviors are as follows (see Appendix A) 



(29) 



with C — 1.303892. The maximal gel size decreases 
as the freezing rate increases. Just below the critical 
freezing rate, the gel is very small as its size shrinks 
quadratically with the distance from the critical point 
5 ~ (1 — q;)^; perturbation analysis shows that this be- 
havior is generic and not limited to the maximal gel size. 
Therefore, freezing is a mechanism for controlling the gel 
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FIG. 2: The maximal gel mass gmax versus the freezing rate 
a (soUd line). Perturbation theory results are shown using 
dashed lines. 



FIG. 4: The multiple gel probability Pmuit versus the freezing 
rate a (solid line). Perturbation theory results are shown 
using a dashed line. 



MASS 




FIG. 3: The mass density m (MASS) versus time r (TIME). 
The system alternates between the coagulation phase and the 
gelation phase. In the former phase the mass decreases lin- 
early according to I20I I such that depletion occurs at time 
r = 1/q. In the latter phase, the active mass decreases slower 
than linear according to 1281 1. The gelation phase ends when 
the gel freezes. 



size: by using freezing rates just below the critical rate, 
it is possible to produce micro-gels. 

The gel freezes following a random Poisson process: 
its lifetime T is distributed according to the exponential 
distribution 



P(T) 



(30) 



As long as the gel is active the system evolves in a deter- 
ministic fashion. When the gel freezes, the total reactive 
mass exhibits a discontinuous downward jump (Fig. 
Since the gel freezes according to a random process, the 
mass of the frozen gel is also a random variable. More- 
over, this quantity is not self- averaging as it fluctuates 
from realization to realization. 

When the gel freezes, the system consists of finite-mass 
clusters only and therefore, the system re-enters the co- 
agulation phase. The system may then undergo a second 
gelation transition that ends when the gel freezes. There- 



fore, the evolution is cyclic with the system alternating 
between the coagulation and gelation. For the same rea- 
son that the gel mass fluctuates, so is the number of 
frozen gels produced a fluctuating quantity. 

The number of gels produced is in principle unlimited, 
i.e., there is a finite probability Pji > that n gels are 
produced. Since the second moment diverges at the gela- 
tion transition according to (|21|l . it is necessarily larger 
than the freezing rate at some finite time interval follow- 
ing the gelation transition. If the gel freezes during this 
time interval then a successive gelation is bound to occur. 
We also note that the evolution in the coagulation phase 
is deterministic and for example, the first two moments 
follow Eqs. (1201 and f^. The "initial" conditions are 
given by the state of the system when the gel freezes. 

It is therefore natural to ask: what is the probability 
that multiple gels are produced? This is the probability 
that the gel freezes prior to time given by Af2(i*) = 
a. Using the second moment M2 = s/(l — sre"*) this 
condition simplifies to 



s(t* 



a e 



-at. 



(31) 



The probability that multiple gels are produced 
is obtained by integrating H30() up to this time, 
-fmuit = Jo* " dT P{T) with the fimiting behaviors (see 
Appendix A) 



Pr 



mult 



0.450851 a T 1- 



(32) 



This probability increases as the freezing rate increases 
(figure 21). It is generally substantial, and moreover, it 
exhibits a discontinuity aX a = a^. 

We now address the final mass distribution of frozen 
gels. Analytic treatment of the successive gelation phases 
is difficult due to the stochastic nature of the freezing pro- 
cess. Numerically, there are two ways to treat the prob- 
lem. One may integrate the rate equations (|28|l up to 
the gel freezing time that is distributed according to H30|) 



7 




FIG. 5: The final mass distribution of the frozen clusters for 
a — 1/2. The simulation results represent an average over 
10^ independent realizations in a system of mass A'^ = 10®. 

and then repeat this procedure if a successive gelation 
occurs. We prefer Monte Carlo simulations where, since 
the system is finite, there is no need to distinguish the 
gel from the finite clusters. In the simulations, we keep 
track of the total aggregation rate Ra = N{Mi — M2)/2 
and of the total freezing rate Rf = aNMo, where N 
is the number of monomers. Aggregation occurs with 
probability Ra/{Ra + Rf)i ^^nd freezing occurs with the 
complementary probability. A cluster is chosen for aggre- 
gation with probability proportional to its mass. Time is 
augmented by At = l/(i?a + Rf) after each aggregation 
or freezing event. 

For the case a < 1, numerical simulations provide con- 
vincing evidence that the tail behavior 

Fkioo) - a < 1 (33) 

is universal (Fig. [SJ. This indicates that frozen clusters 
produced during the coagulation phase dominate at large 
masses as every such phase is expected to contribute 
to the tail according to Eq. (|27ll . Intuitively, this is clear 
because large clusters are quickly merged into the gel and 
therefore, frozen clusters produced in the presence of the 
gel tend to be small. Interestingly, freezing leads to a 
new critical exponent in mean-field percolation. 



V. CONCLUSIONS 

In summary, we studied aggregation processes with 
freezing. For constant freezing rates, the problem can 
be formally reduced to the no freezing case. The mass 
distribution of frozen clusters resembles the mass dis- 
tribution of reactive clusters, decaying exponentially at 
large masses, when the freezing is sufficiently fast. Novel 
behavior emerges when the freezing rate is slower than 
some critical value. In this case, the mass distribution of 
frozen clusters decays algebraically. For constant aggre- 
gation rate, the critical freezing rate is zero but for the 
product aggregation rate, it is finite. 

For the product aggregation rate, the freezing rate con- 
trols the gelation process. If it is sufficiently high, no 
gelation occurs, and if it is just below the threshold, 
micro-gels are produced. If one gel is produced, then 
multiple gels are possible. In this case, the mass of the 
gels and their number are both controlled by a random 
process and as a result, they fluctuate from realization 
to realization. The system exhibits a series of gelation 
transitions and it alternates between ordinary coagula- 
tion and gelation. The random freezing process governs 
the number of percolation transitions as well as the mass 
of the frozen gels. 

The behavior when freezing occurs spontaneously is 
quite different than that found when freezing is reaction- 
induced (upon merger a cluster may freeze with some 
fixed probability) [2^ l^^. For the constant aggrega- 
tion rate, the mass distribution is always algebraic and 
the characteristic exponent is non-universal as it depends 
on the freezing probability. For the product aggregation 
rate, the gelation transition is always suppressed, because 
the probability that a cluster remains reactive decays ex- 
ponentially with the number of merger events. 

There are many open questions raised by this study. 
For example, it will be interesting to investigate the be- 
havior in low-dimensional systems where the rate equa- 
tion approach no longer holds. Additionally, the expo- 
nent characterizing the mass distribution of frozen clus- 
ters represents a novel critical exponent in percolation 
processes and this should be a challenging problem be- 
low the critical dimension. 
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APPENDIX A: PERTURBATION ANALYSIS 
1. a T 1 

To investigate the behavior in the time domain 1 < r < 
a~ ^ , we make the transformation t = 1 + (q:^^ — l)x with 
< X < 1. The governing equations (|28|) become 



dm 
dx 
ds 
dx 



1-x 
s 

1-x 



s{m — s) 



ail — X 



(Ala) 
(Alb) 



sx 



The initial conditions are now m(0) = s(0) = 1 — a. 
When at 1, we perform perturbation analysis using the 
small parameter t = 1 — a. Given the initial conditions 
to(0) = s(0) ~ e, we write 



s = eGi 



(A2) 



Substituting these expansions into the governing equa- 
tions (jAla|) - (jAlb|) and keeping only dominant, linear in 
e, terms we obtain 



dFi 
dx 
dGi 
dx 



1-x' 
Gi 
1-x' 



(A3a) 
(A3b) 



The initial conditions are Fi{0) = Gi{0) — 1. Solving 
these equations, we find 



Fi^Gi 



1 - x. 



(A4) 



Thus, the gel mass vanishes to first order in e, and we 
should consider the second order terms. The second order 
terms are coupled according to 

dFj. _ G2 

dx 1 — X 

F2 — G2 



dG2 ^ G2 _|_ n _ 
dx ^ 1-x ^ G2+ x{l - x)' 



(A5a) 
(A5b) 



The boundary conditions are F2(0) = G2(0) = 0. We 
seek the non-trivial solution with the following deriva- 
tives at the origin: F^(0) = 0, and £"2(0) = -2. 

To find the probability that multiple gels are produced, 
we notice that the Poisson distribution 1)30(1 is uniform 
in terms of the transformed time, P{x) — 1 and there- 
fore, the probability of forming multiple gels is simply 
^muit = Jq* dxP{x) — x*. Also, the condition be- 
comes s(x*) = a(l — a)(l — a;,). Substituting the pertur- 
bative expansion e(l — x*) -I- e^G2(x,) — e(l — e)(l — x,), 
this condition becomes G2(x*) = x, — 1. It also implies 
62(1) = 0. Numerical integration yields -Pmuit = x» = 
0.45081. 

To find the maximal gel mass, we write g — m — s ^ 
£^(^2 — G2) and use 62(1) = 0. Thus, the maximal gel 
mass g,nax = g(x = 1) is 



9n 



(1 



(A6) 



Integrating these equations numerically, we determine 
^2(1) = 1.303892. Moreover, the perturbation analysis 
shows that in general, the gel mass vanishes quadratically 
close to the critical freezing rate, g oc (1 — a)^. 



2. a [0 

When freezing is slow, we may drop the freezing loss 
term —aRk from the master equation H19|) as done in sec- 
tion III. The problem therefore reduces to the no-freezing 
case where u = t = t and thus, the generating function 
^ is 3l(z,i) = t-^G{z + Int - t). Invoking the iden- 
tity Ge^'^ = e/- , the sol mass s = 3l(z = 0) satisfies 



In the long time limit, 
s{t) = e"* + te"2t ^ 



(AT) 



From the one-gel criterion 

and (HZl) we obtain t.^, — 
In i. Finally, using P^^yt = 1~ exp[-a(t, - 1)], we 
find the multiple-gel probability: Pmuit — — Q;ln(ea). 

To leading order, the total mass m remains constant. 
Using the exact governing equation dm/dt = —as, we 
derive the first order correction: 



m(t) = 1 - a - a / dt' s{t'). 



The maximal mass of the gel is therefore 



5max = 1 - Ba, B =1 



dts{t). 



(A8) 



(A9) 



Using t = —{1 — s) ^ In s we change the integration vari- 
able and then transform the integral: 



B = 1 



ds s 



1 



In; 



S{l-S) (1-S)2 



(AlO) 



Performing the integration, we find B — 7r^/6. 



